Birkhoff normal forms for Fourier integral operators II. 



A. lantchenko* and J. Sjostrandf 

Abstract. In this work we construct logarithms and Birkhoff normal forms for elliptic Fourier integral operators 
in the semi-classical limit, under more general assumptions than in a previous work by the first author. The 
methods are similar but slightly different. 

Resume. Dans ce travail nous construisons des logarithmes et des formes normales de Birkhoff pour des 
operateurs integraux de Fourier dans la limite semi-classique. Les hypotheses sont plus faibles que dans un 
travail anterieur du premier auteur et les methodes sont semblables mais un peu differentes. 

0. Introduction. 

This work is a continuation of the work [la] of the first author. As in that paper we 
consider the problem of constructing a microlocal logarithm and a normal form of an eUiptic 
semi- classical Fourier integral operator near a fixed point of the corresponding canonical 
transformation. In [la] the canonical transformation was assumed to be of real hyperbolic 
type and the purpose of the present work is to relax this assumption considerably, to what 
we think are the natural conditions. 

As in [la] a motivation is to improve the analysis of quantized billiard ball maps near 
closed trajectories for boundary value problems. Then the associated canonical transfor- 
mation is the corresponding classical Poincare map. In the introduction of [la] one such 
problem in the case of scattering by obstacles was mentioned, where this canonical trans- 
formation is of hyperbolic type near its fixed point, but in many other cases the canonical 
transformationcan be more arbitrary. (Such an improvement of spectral asymptotics for 
non-degenerate potential wells was obtained in [Sj] by means of a Birkhoff normal form.) 

Another motivation for returning to this study comes from recent works on inverse 
spectral problems by Guillemin [Gui] and Zelditch [Zel-3]. In these works, quantum 
normal forms of the whole (Laplace) operator were constructed in a neighborhood of the 
closed trajectory. We believe that it may often be sufficient and technically easier to work 
with normal forms of a corresponding monodromy operator, which coincides with the 
quantized billiard ball map in the case of boundary problems. For general semi-classical 
problems, this operator was recently introduced and studied in a more systematic fashion 
in [SjZw]. In particular one might arrive at very nice trace formulae by combining the 
basic trace formula of that paper with the normal form obtained in the present work. 

It turned out to be somewhat difficult to extend the whole method of [la] to the 
more general situation here. In [la] a major idea was to use that the symplectic group is 
connected and work with deformations from the the identity transformation to the given 
canonical transformation. In the hyperbolic case this can be done in a such a way that the 
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intermediate transforms satisfy the assumptions for having a normal form. In the general 
case considered here, we encounter singular values for the deformation parameter where the 
conditions are not fulfilled. It seems possible to circumvent this difficulty by complexifying 
the deformation parameter and use corresponding slightly complex transforms. Eventually 
however we found a method allowing us to avoid deforming the differential of the canonical 
transformation but only the higher order part in its Taylor expansion. In this way we have 
families of objects which satisfy our assumptions everywhere. 

Consider a semiclassical Fourier integral operator A of order 0, with an associated 
canonical transformation k : neigh ((0, 0), R^"^) — > neigh ((0, 0), R^"^) having (0,0) as a 
fixed point. Assume A is elliptic at (0,0). The set of eigen- values of dK.{0) is then closed 
under inversion A — 1/A and under complex conjugation. Assume that to the distinct 
A in the spectrum of dK(0,0), we can associate a logartithm ^ = log A in such a way 
that inversion and complex conjugation correspond to the map fi — |U and to complex 
conjugation respectively. (Notice that this assumption excludes the existence of negative 
eigen- values of dK,{0,0).) Assume also that 

kjfij = 27rm, kj, n e Z =^ kjfXj = (0-1) 

Then a real version the Lewis-Sternberg theorem (see [St], [Fr] and Theorem 1.3 below) 
tells us that we can write k{p) = expHp{p) + 0{p°°) {p = {x,^) e R^") for some smooth 
and real- valued p = 0{p'^). The first result of this work says that under the same as- 
sumptions, we can write A = e~*^/^ modulo an operator which vanishes to infinite order 
at (x,^) = (0,0), h — 0, where /i > is the small semi- classical parameter (Theorem 
3.2). Here P is a semi-classical pseudodifferential operator of order with p as its leading 
symbol. 

The second result is a straight forward extension of the normal form in [Sj] and says 
that under the non-resonance condition (4.2) below, the operator P has a simple normal 
form in terms of certain action operators. (In (0.1) we do not have to be very precise 
concerning the enumeration of the logaritms of the eigen-values of the linearization as long 
as we have have one of each, modulo the sign, while in (4.2) we enumerate one rather 
specific half of these logs. Also notice that the combination of the two conditions takes 
the simple form (4.3).) We do not expect that the exclusion of negative eigen-values is a 
serious restriction, for if such eigen-values are present, we can find specially adapted and 
explicit metaplectic operators M and apply our results to MA. 

Both our results were obtained by the first author ([la]) under the assumption that 
dK{0, 0) is of real- hyperbolic type. 

For the reader's convenience we have taken the pain review some well-known linear 
symplectic geometry in section 1, and in that section we also give a (probably well-known) 
proof of the real version of the Lewis-Sternberg theorem, which has a structure that we can 
follow in the proof of the corresponding quantum result for logarithms of Fourier integral 
operators. 

It is beyond the scope of this work to consider convergence questions related to the 
perturbation series that appear in connection with normal forms. See for instance H. Ito 
[It] and G. Popov [Pol, 2]. 

The plan of the paper is the following: 
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In section 1 we review some standard facts about linear symplectic geometry and add a 
few remarks for treating the real case. We also review a proof of a real version of the 
Lewis-Sternberg theorem, that we can later use as guideline for the proof of the quantum 
result. 

In section 2, we introduce some notions of equivalence that are used in the formulation of 
the main results. These notions are essentially the same as in [la]. 

In section 3 we establish the main result about logarithms of Fourier integral operators. 

In section 4 we give the "Birkhoff" normal form for the logarithm i.e. for a certain h- 
pseudodifferential operator of order 0. 

In section 5 we extend the results to the parameter dependent case. In many genuinely 

semi-classical problems we do not have any homogeneity, inferring that the Poincarc map 
is essentially energy independent, and then we cannot expect in general that our arithmetic 
condions be fulfilled for all energies in some interval. Consequently it is of interest to know 
that the results are valid to infinite order at points where the conditions are fulfilled. 

1. Review of some symplectic geometry. 

In this section, we review some more or less well-known arguments that will later be 
extended to the quantized case. See [MeHa]. 

We start with the linear case and let A : C^" —>■ C^" be linear and symplectic in the 
sense that 

a{Ax, Ay) = a{x, y), x,y e C'", (1.1) 

where a is the standard symplectic 2-form on C^"^. We recall that this implies that 
det^ = 1, since A will conserve the volume form a'^/nl. When n = 1 the converse is 
also true. 

Let Ex = AfdA — A)^") = Ker ((^ — A)^"^) be the generalized eigen-space associated 
to A e C \ {0}. 

Proposition 1.1. If Xfi ^ 1, X, ji & C\{0}, then Ex and E^ are symplectically orthogonal: 
Ex ±" E^. 

Proof. Let E^^^ = ExnAf{{A- A)-?), so that 

^ E^^^ c E^^^ c .. C E?") = Ex. 

Define Ejj^'^ in the same way. Then for x e E^\ y e Ej^^: 

111 1 

a{x, y) = a{-Ax, -Ay) = —a{Ax, Ay) = j^cr{x, y), 

and since l/X/j, 1, we get a{x, y) — 0. 
Assume that we have for some m > 2: 

a{x,y) = 0, for x e eI^\ y e E^^^ j + k < m. (P^) 
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We have just established (P2)- 

Let X e E^^^ , y e Ef \ j + k = m+ 1. Write A\^^ = X + N, = ji + M, where 

N,M are nilpotent with N : e[^^ e[^~^\ M : E^j"^ E\jt~^\ Then 
a{x, y) = ai^Ax - \nx, -Ay - -My) 

A A fJ, fl 

= -^aix, y) + (j{-^x, -Ay) + a{\Ax, --My) + a{\Nx, -My). 
AfJi A jJi A jJi A fJj 

The last three terms vanish because of the induction assumption, and we get a{x,y) = 0, 
so we have proved (P^_|_i). The proposition follows. # 

We conclude that Ex are isotropic if 7^ 1 (i.e. cr vanishes on Ex x Ex). We also see 

that 

El ±^ Ex for A 7^ 1, E-i E^ for ^ -1. 

It follows that 

k 

= El e E_i e 0(£;a, e (1.2) 

1 ^ 

where A^, A~^ and possibly 1,-1 denote the distinct eigen-values of A with Aj 7^ ±1. 
Moreover, all the © and except the ©s inside the parenthesies indicate symplectically 
orthogonal decomposition. 

For A ^ {1, — 1}, write A\^^ = A + Nx with Nx nilpotent. a\^^^^^ is non-degenerate 

and if '^A denotes the symplectic transpose of A, so that a{Ax.,y) — a{x,^Ay), then 
'^A = A, and "AA = 1, and hence 

''A = A-\ 

Also notice that the Ex are invariant subspaces for "'A. On Ei/x we have on the one hand 
A= 1/X + Ni/x and on the other hand, 

1 = -4^ = (A + -iVA)(i + iVi ) = (1 + j''Nx){l + AiVi ). 



Hence, 



AiVi = -ij'^Nx) + (^^TVa)' - (..)' + - • (1-3) 



On El we have A = 1 + A^i, where A^i is nilpotent. The requirement that A be 
symplectic means that 

(l + '^iVi)(l + iVi) = l. (1.4) 

Similarly on E_i, we have 

C^iV_i-l)(iV_i-l) = l. (1.5) 

Conversely consider a decomposition of C^" as in (1.2) into symplectically orthogonal 
spaces El, E^i, Ex^ ®Ei/Xj with E±i symplectic, Ex^, -E'l/Aj isotropic for A-,- 7^ 1, —1 and 



all the 1, — 1, Aj,l/Aj different. Let A be an operator leaving all the £■( ) invariant, with 
= A + Nx, A = ±1, Aj, 1/Xj and N\ nilpotent. Then A will be symplectic if (1.3-5) 

hold. 

We next consider logarithms of a symplectic matrix A. Decompose C^" into general- 
ized eigen-spaces as in (1.2). We construct log^ with the same generalized eigen-spaces 
in the following way: On Ei we have A = 1 + Ni with A^i nilpotent and we put 

log^ = log(l + iVi) = iVi - ^Nf + ^Ni + 

where the sum is finite, and in the following we always define the log of 1 + in this way, 
when N is nilpotent. If A G {Aj, A~^, —1}, we choose /x = A*(A) with A = e'^ in such a way 
that 

f^ih = -/^(A,). (1.6) 

Write 

A\^^=X + Nx = X{l + jNx), 

and define on Ey. 

\ogA = ii + \ogil + jNx). 

This gives a definition which only depends on a choice of logarithms of Xj, 1 < j < k, and 
on log(— 1) if E_i has positive dimension. It is easy to check that 

log^ + '"log^ = {2k + l)27rzn_i, (1.7) 

for some integer k, where n_i denotes the spectral projection onto E-i. Of course we 
have that exp log A = A. Recall also that ii B + "^B = 0, then exp i? is a symplectic matrix. 

Assume now in addition that ^ is a real matrix: A : R^" — > R^*^. Then E±i become 
real in the sense that they are invariant under complex conjugation T : (x,^) i-^ {x,^). 
The same holds for Ex if A is real. If A is not real, there are two possibilities: 

1) |A| 1. Then A, j, = are also eigen- values and 

Ex® El ® Ey® El 

is the complexification of a real symplectic space. 

2) |A| = 1. Then 1/A = A and Ex ® Ei is the complexification of a real symplectic space. 

In all cases we have 

and it is easy to see that log^ will enjoy the same property, provided that we choose /i(A) 
in such a way that 

/.(A) = MA). (1.8) 
This is possible, if we assume that A has no negative eigen-values. We get 
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Proposition 1.2. a) Let A be a complex symplectic 2n-matrix, and choose a value //(A) 
of the logarithm of each distinct eigen-value X, different from 1 in such a way that (1-6) 
holds. Then we have a corresponding choice of B = log A with e^ = A, satisfying (1.7). 

b) Assume in addition that A is real and has no negative eigen-values. Then by choosing 
n{X) with the additional property (1.8), B = log A becomes a real matrix, and '^B + B = 0. 

From now on, we work under the assumptions of b) above, so that B = log A is real 
and symplectically anti-symmetric. Consider the quadratic form 

b{p) = ^a{p,Bp). (1.9) 

Then 

{db{p),t) = ^{a{t, Bp) + a{p, Bt)) = a{t, Bp). 

On the other hand 

{db{p),t) = a{t,Hi,{p)), 
where Hi) denotes the Hamilton vector field associated to the function b, so 

Hbip) = Bp, (1.10) 

and the fact that B = log A can be expressed by 

A = expi/b. (1.11) 

We now consider the problem of finding the " logarithm" of a canonical transformation 
also in the non-linear case. We first proceed somewhat formally and let Ps be a smooth 
real function depending smoothly on the real parameter s. Consider the corresponding 
canonical transformation 

Kt,s = (iWtHp^. (1-12) 

We will later assume that ps vanishes to second order at some point po, and then the germ 
of Kt,s at Po will be well-defined for all real t. We differentiate the identity 

dtKt,s{p) = Hp^{Kt,s{p)), 

with respect to s: 

dtid.ntM) - {-^{KtM){d,ntM) = {dsHpJiKtM) = Ho^p^iKtM). (1.13) 

Notice that the differential dKt^s{p)i^ — ^gp" (p)^ satisfies 

dH 

dtdKt,s{p) Q^{Kt,s{p)) o dKt,s{p) = 0, dAvo,s(p) = 1. (1.14) 
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Comparing the last two identities, we see that 

ft 



dsKtM = f d{Kt_tJ{Ki^sip))iHd,Ps)i'^iAp))dt, (1.15) 

^0 



which can also be written as 

ft 



daKt,a = f )^{Hd^pJdt, (1.16) 

^0 



where we use standard notation: lower * for push forward and upper * for pull back. 
Rewrite (1.16) as 



Jo 



and notice that 

since ^ is a canonical transformation. Then 



Jo ' Jo 



SO 

dsKt,s ^ {Kt,s)*Hq^^, (1.17) 

where ^ 

Qt,s= dsPsOKfgdt. (1.18) 
Jo 

In the last formula we shall take t — 1 and consider a problem where dgPs will be the 
unknown. More precisely, let k be a smooth canonical transformation: neigh (0, R^") — > 
neigh (0, R^"^) with k{0) = 0. Let A := dK,{0) =: kq have no negative cigen-values so that 
part b) of Proposition 1.2 applies. (More assumptions will be added later.) Let i? be a 
real logarithm of A as in the proposition and define the quadratic form po = b as in (1.9). 
Then 

Ko = expHp^. (1-19) 
We look for p{p) — po{p) + 0{p^), so that 

«(p) = expHpip) + 0(p°°). (1.20) 

Let «;s, < s < 1, be a smooth family of canonical transformations with 

Ks{0) = 0, dKs{0) = cZk(O), ko = cZk(O), ki = k. (1.21) 

Then we look for a corresponding smooth family Ps{p) — Po{p) + C'(p^), with Ps=o — Po as 
above, such that 

Ks{p) = expHpSp) + 0{p°^). (1.22) 



Then p = pi will be a solution to our problem. Define qs{p) = 0{p^), by dsKs = {Ks)*Hq^, 



or: 



{KsydsKs = Hg^. (1.23) 



The discussion leading to (1.18) indicates that we should find ps with the above properties, 
so that 



i(p) = / dsPs 
Jo 



oexptHp^{p)dt + 0{p°°). (1.24) 



Let N > 2 and suppose that we have already found a smooth family pi^\p) = 
Po{p) + 0{p^) with =po, so that 

qs{p)= [ a,p(^)oexpti? w(p)dt + e»(p^+i). (En) 
Jo 

Notice that p^^^ = po solves (E2) since qs{p) = 0{p^). Look for pi^'^^'* = pi^^ + ri^'^^\ 
with rf+^\p) = 0{sp^+^). Then 

exptH (N+i) (p) = exptH (N) (p) + 0{p^), 
e^ptH (N+i){p) = exptHp^ip) + O(p^), 



and we get 

■ 1 

a,p(^+i)oexptiy (^+i)(p)(it = 

C o exptH <^,{p)dt + C dsTf+^^ o exptHp,{p)dt + 0{p^+^). 

Jo Jo 



f 

Jo 



If we write the remainder in (Eat) as —v^^''^^\p)+0{p^~^'^), where is a homogeneous 

polynomial of degree + 1 depending smoothly on s, we will get (Ejv_|_i), if we can find 

dgri^'^^'^ as a homogeneous polynomial w*^^"*"^-* of degree A'' + 1 depending smoothly on s, 
such that 



i 



1 



o exptHp^{p)dt = v^^+^\p). (1.25) 



'0 

Consider a general linear map B : C^" — > C^" with Jordan decomposition B = D + N , 
where D is diagonalizable, = diag {dj) with respect to a suitable basis, and N is nilpotent 
and commutes with D. The action of the vector field Bx ■ dx on the space (C^"^)* of linear 
forms on C^"^ can then be identified with *B in the natural way. Notice that *B has the 
Jordan decomposition *D +*A^ and that *D becomes diag {dj) if we select the dual basis to 
the one where D is diagonal. Define 

tg(m) = ®l® ..®l + l®*B ®l® + ...1 (8) .. (g) 1 *B, 
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as a linear endomorphism of the m-fold tensor product ((C^")*)®"*. We have the Jordan 
decomposition 

(where the first term to the right is diagonahzable, the second nilpotent and the two 
terms commute). The corresponding eigen- values are dj^ + .. + dj^, for j = [ji, ..,jm) £ 
{1,2, ..,2n}^^''^'"''^^. The three operators in (1.26) act naturally on the symmetric tensor 
product (^(^c2n^*^0m ^y^^ decomposition (1.26) is still a Jordan one on that space. The 
eigen- values of ^B^'^^ become 

Pk = kidi + .. + k2nd2n, k e N^'^, ki + .. + k2n = m. (1.27) 

^^Q2ny y3m Qq^oi to the spacc 7^hom(^^") TO-homogcneous polynomials on C^" and 
is the action of Bx ■ on that space. 
Consider the map 

^0 

which is equal to 

/ eMt'B^'^^)dt. (1.29) 
Jo 

The Jordan decomposition (1.26) gives a similar decomposition of (1.29). The eigen- values 
of (1.29) are therefore given by 

1 r 1 for = 

6'^"^= { e^'^ -1 ^ ^ , ^ , /c e N^", l/cl =m. (1.30) 

^ — ^ for /3fe ^ ' 

Pk 







We conclude that the map (1.28) is invertible for a given m precisely when for all k G N^" 
with \k\ = m: 

kidi + .. + k2nd2n £ 27rzZ =^ kidi + .. + k2nd2n = 0. (1-31) 

Now return to the equation (1.25), where po{p) = b{p) = ^a{p,Bp) and B is the 
logarithm of the real symplectic matrix A = dK{0), obtained under the assumptions of 
Proposition 1.2, part b). The eigen- values of B are then of the form with some possibly 
vanishing multiplicity and ^ with equal multiplicity > 0. Here we arrange so 

that all the eigen- values are distinct for instance by taking //^ with either Re pj > or 
with Re pj = and Im pj > 0. We also recall that our set of eigen- values is closed under 
complex conjugation. The assumption that (1.31) holds for all m, then amounts to the 
assumption that 

^ kjPj e 27riZ =^ ^ kjPj = 0, (1.32) 

for all ki, ..,kr G Z. Here r < n is the number of distinct pj. (We could also have chosen 
to repeat the eigen-values according to their multiplicity without changing (1.32).) 
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We have practically finished the proof of the following version of a theorem of Lewis- 
Sternberg: 

Theorem 1.3. Let k : neigh (0, R^"') neigh (0, R^") be a smooth canonical transfor- 
mation. Assume that dK{0) has no negative eigen-values. Let the distinct eigen-values of 
dK{0) be 1 (possibly with multiplicity 0) and Xj, A~^, 1 < J < r with \Xj\ > 1 or with 
\Xj \ — 1 and < argAj < n. Choose nj with Xj = e^^ , in such a way that Xj corresponds 
to Jtj and let B = log A be given as in part b of Proposition 1.2. Assume that (1.32) holds, 
and let po{p) = b{p) be given in (1.9). 

Then there exists p{p) e C°° (neigh (0, R^"); R) such that p{p) = po{p) + 0{p^) and 

K{p)^expHp{p) + 0{p'^). (1.33) 

p is uniquely determined by these properties (for a given choice of po). 

This result (at least the existence part) is extremely close to a corresponding one for 
complex canonical transformations, due to Lewis-Sternberg ([St], Theorem 1 and Corollary 
1.1) and clearly stated in [Fr], Theoreme V.l. 

End of the proof. We establish the existence of p. Let Kg, < s < 1 be a smooth 
family of canonical transformations with (0) = 0, dKs{0) = dn{0) and with kq 

linear (= dK,{0)). Define qg by (1.23). The preceding discussion gives us a smooth family 
Ps{p) e C°° (neigh (0, R^"); R) with ps=o = POi such that if Kg = expii/p^, then 

K:dgKg^H,^+0{p°^). (1.34) 

A simple computation shows that (1.23,34) can be written as 

dgKj^p) = -H,^{k;\p)), dgK-\p) = -H^SK\p)) + O(p-), (1.35) 

and we also know that Kq^ =Kq^. It follows that 

K-\p)^n-^ + 0{p'^), 

and hence that 

Ks{p) = Ks{p) + 0{p'^). 

Taking s = 1 gives (1.33) with p = pi. 

We next prove the uniqueness of the Taylor expansion of p. Let p have the same 
properties as p, so that 

exp HM = exp H^p) + 0{p^), p{p) = p{p) + 0{p^). 

Assume that p — p does not vanish to infinite order, so that p = p+r, where r(p) = 0{p^), 
r{p) ^ (9(p"^+^), for some 3 < m G N. 

Put Ps — {1 — s)p + sp — p + sr, < s < 1, so that po — p, Pi — p, and define 
Kg = exp Hp^ , so that 

«i(p) = «o(p) + 0(p~). (1.36) 
10 



For this family, define Qs by (1.23). Then (1.24) holds and since dsPs = r, we have 

qs= f roexp(tiypJdt + 0(p^+i). 

The previous discussion shows that the integral has a non-zero Taylor polynomial of degree 
m: 

qsip) = q{p) + 0(P"^+'), ^ g e prom- (1-37) 
From (1.35), we conclude that 

which contradicts (1.36). The proof is complete. # 

2. Notions of equivalence. 

As in [la], our results will be valid "to infinite order at (0,0)" and in this section we 
review the corresponding notions of equivalence. Using these notions we also develop a 
very rudimentary functional calculus for functions of several pseudo differential operators. 

If Vj C R"^ are open neighborhoods of and vj G C°°(y^), we say that vi and V2 
are equivalent; vi = V2, if vi — V2 vanishes to infinite order at 0: vi{x) — V2{x) = 0{x°°). 
This is clearly an equivalence relation and the equivalence classes can be identified with 
the corresponding formal Taylor expansions. 

With V = Vj G R"^ as above, let SI^y{V) denote the space of functions a{x]h) in 
C°°{V) depending on the semi-classical parameter h g]0, ho] for some ho > 0, such that 

oo 

a{x,h) ^^h^aj{x), h^O, (2.1) 



for some sequence of aj G C^{V). We say that a^'^^ G S'^j(Vfc), A; = 1, 2 are equivalent and 
write a*-^-* = a^'^\ if a^-^"* = a^-^'* for the corresponding coefficients in (2.1). Equivalently we 
can say that a^'^^ are equivalent if a^^\x; h) — a'^'^\x; h) — 0{{x, h)°^). 

If m > is a smooth weight function on V, we define S'^(V, m) to be the space of 
smooth functions a on F such that for all multi- indices a, we have \d^a{x)\ < Ca,a'm'ix)- 
Let 5'°i(y,m) be the subspace of S^^ for which (2.1) holds in S^{V,m). 

We next pass to the case of pseudodifferential operators. Recall that if p{x, ^) belongs 
to an appropriate symbol class of functions on R^"^, then we define the corresponding 
/i-Weyl quantization P = p^{x, hD^) by: 

^""^"^^ = (2^ / i e)u{y)dyde. (2.2) 

Recall that p is called the Weyl-symbol of P. (See for instance [DiSj].) Let S'°(R^'^) denote 
the space of smooth functions that are bounded together with all their derivatives. If 
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p^^'^ e (5°, n 5°)(R2^) = ^.V^-^"' 1), ^ = 1, 2, we say that (p('=))"'(a;, hD; h) are equivalent 
(and use the symbol =) if p^^^ are equivalent in the sense of the classes S^^. 

We will use the abbreviation neigh (0, R") to denote some neighborhood of in R". 
Tacitly it is understood that these and other geometrical objects are independent of h. We 
say that two smooth canonical transformations Kj : neigh (0, R^"') — > neigh (0, R^"^) with 
Kj{0) = are equivalent if Ki(p) — K2(p) = 0{p°°). Possibly after shrinking the neighbor- 
hoods, we can introduce the inverses kJ^ : neigh (0, R^"^) — > neigh (0, R^"^). Then ki = K2 
iff = i^'2^ ■ Also notice that the notion of equivalence of canonical transformations is 
stable under composition in the natural way. 

Let K : neigh (0, R^"^) — > neigh (0, R^") be a canonical transformation which maps 
to 0. Then there exist N e 'N and a non-degenerate phase function 0(a;, y, 9) e 
neigh (0, R"^+'^+^) such that the graph of k in a neighborhood of (0, 0) coincides with 
the image of the local diffeomorphism: 

:= {(x, y, 0); c^q{x, y, ^) = 0} 9 (x, y, d) ^ (x, c/)',; y, -</.;). (2.3) 

Here we recall that a smooth real-valued function is called a non-degenerate phase function 
(in the sense of H5rmander) if d(f)'g^, .., d(f)0^ are linearly independent on the set C(j, above, 
which then becomes a 2n-dimensional smooth sub-manifold. When discussing the relation 
between phases and symbols with canonical transformations, it is tacitly understood that 
the point a; = 0,y = 0,^ = corresponds to k(0) = under the map (2.3). 

Let K be as above and let be a corresponding generating phase. Let k : 
neigh (0, R^"^) neigh (0, R^"^) be a second canonical transformation (with the tacit con- 
vention that it also maps to 0). It is easy to see that k = k if and only if k has a 
generating phase (f) which is equivalent to (f). 

With 0, K, as above we consider a Fourier integral operator of order 0: 

Uui.) = 0; KHvHy^, (2.4) 

where a G 5*^1 has its support in a sufficiently small neighborhood of (0, 0, 0). In this paper 
we only consider Fourier integral operators that are elliptic at (0, 0, 0) in the sense that 
ao(0, 0, 0) ^ 0. In order to normalize things, we will always assume that 0(0, 0, 0) = 0. If 
K is the canonical transformation generated by 0, we say that k is the canonical transfor- 
mation associated to U . Thanks to the ellipticity assumption, k, is uniquely determined 
by U in some neighborhood of 0. We also recall the fundamental theorem about Fourier 
integral operators, namely that if '^(a;, j/, w) is a second phase which generates n and if a 
has support in a sufficiently small neighborhood of (0,0,0), then there exists a classical 
symbol b{x,y,w;h) of order with support in a small neighborhood of (0,0,0), such that 
/(6, '0) (formed as in (2.4) with N replaced by the dimension of w-space) is equal to I (a, 0). 
Let K be a second canonical transformation with 'k = k. Let (p = (p he a, corresponding 
generating phase. We say that U = I{a, (p) is equivalent to U and write U = U,ifa = a. 
It is a standard exercise in Fourier integral operator theory to verify that this definition of 
equivalence does not depend on the choice 0. It is also easy to show that the definition is 
stable under composition in the natural way. 
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Below we will also need some functional calculus. First we consider exponentials 
of pseudodiffcrential operators. Let p ~ Yl^ Pji^^O^'' 5'^(R^"^,1), and assume that 
Po is real-valued with po(0,0) = 0, j3o(0,0) = 0. Then e'^^P"^'^''^^''^^/'^ is well-defined 
for all complex t (even without the reality assumption on po) and for real t we get a 
Fourier integral operator. If x £ C^(R^") is equal to 1 near 0, then up to an operator 
whose distribution kernel is rapidly decreasing together with all its derivatives, we have 
that ■^'^Q-'^tP/h^w jg ^ Fourier integral operator as above, with the associated canonical 
transformation nt = ex^ptHp, whose equivalence class does not depend on the choice of x- 
It is also easy to see that if P is a second pseudodiffcrential operator which is equivalent 
to P and with real leading symbol, then for real we have e~**^/'* = g"**^/'* (in the sense 
that we have equivalence for the corresponding truncated operators). 

Finally we discuss a very primitive pseudodifferential functional calculus. Let Pk = 
Pk{x, hD; h), k = 1, .., A'^o be a commuting family of pseudodifferential operators with pk G 
S'°j(R2'^, ((x,0)"') (with the standard notation ((x, ^)) = (l + |(x, V2) and assume that 
the leading symbols pk^o vanish at (0, 0). Let F{li, .., iNo\h) e 5'^j(neigh (0, R^")). Let Fjv 
be the sequence of polynomials in li, .., ln^, h obtained by taking the Taylor polynomials 
of order N of the first N terms in the asymptotic expansion of F, so that 

F-FAr = 0((6,/i)^), {L,h)^0. 

Then it is easy to see that Fn{Pi, ••, Pn, h) = q^{x^ hDx'-, h), where 

QNix, t h) - qM{x, t h) = 0{{x, hf^^^^^), where k{N, M) oo, iV, M oo, 

and that this sequence defines naturally an equivalence class of pseudodifferential operators 
that we shall denote by F{Pi, .., Pnq'i h)- 

3. Logarithms of Fourier integral operators. 

Let ?7s, < s < 1 be a smooth family of elliptic Fourier integral operators of order 
0, associated to a fixed canonical transformation k : neigh (0, R^") — > neigh (0, R^"') with 
k(0) = 0. We represent Ug by 

U,u{x) ^ JJ e^'l'^'''y^'^^Us{x,y,e;h)u{y)dyd9, (3.1) 

where Ug G 5*^1 and more generally d^Us G S^^ for all /c G N is a smooth family of classical 
symbols of order 0, defined in neigh ((0, 0, 0); R^"'+-'^) and </> is a real phase function which 
is non-degenerate in the sense of Hormander [Ho] and generates k, so that C<^ 3 {x,y,9) ^ 
{x, (f)'^; y, —4>'y) ^ graphs is a local diffeomorphism, where Crf, C R^n+Ar -g ^j^^ sub-manifold 
given by 0^(a;, j/, 9) — 0. To normalize things, we assume that 0'(O, 0, 0) = and that 

0(O,O,O) = O. (3.2) 

Notice that this last assumption does not depend on the choice of phase in the repre- 
sentation (3.1). In the following, we shall use the equivalence relations "=", defined in 
section2. 
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We define the "logarithmic derivative" of our family, to be the smooth family of 
pseudodifferential operators Qs given by 

Qs = U-^hDsU,. (3.3) 

Qs and more generally d^Qs is a smooth family of classical pseudodifferential operators 
defined in neigh ((0, 0); R^"). (We made an arbitrary choice of the order of the factors in 
(3.3), if Qs = {hDsUs)U~^ , then we get a new pseudodifferential operator which is related 
to Qs by the intertwining relation UsQs = QsUs-) 

The family Us is determined uniquely by Uq and its logarithmic derivative: 

Lemma 3.1. Let Vg be a second family of Fourier integral operators with the same 
properties as Us and associated to the same canonical transformation k. Assume that 
U-^hDsUs = Vs-^hDsVs and that Uq = Vq. Then Us = Vs. 

Proof. Let U he a fixed elliptic Fourier integral operator associated to k, so that 

Us = UAs, Vs = UBs, 
where As, Bs are smooth families of pseudodifferential operators. Then 

U-^hDsUs = Aj^hDsAs, 

and similarly for Vs, so we get 

A;^hDsAs = B-^hDsBs, Ao = Bo. (3.4) 

^Prom this we conclude first that As and Bs have equivalent principal symbols, then 
equivalent sub-principal symbols and so on, so As = Bs and hence Us = Vs. # 

Remark. We have HDsUs = UsQs, hence HDsU* = —Q^U* for the adjoint operators, so 

hDs{u:us) + iQ:{u:us) - iu:us)Qs) = o. 

UgUs is a smooth family of elliptic pseudodifferential operators and we conclude 

a) If Us are unitary (up to equivalence), then Qs are self-adjoint (up to equivalence). 

b) If Us is unitary for one value of s and Qs are self-adjoint for all s, then Us is unitary 
for all s (again up to equivalence). 

We now assume for a while that 

Us = Ui,s, where Ut,s = e-**^«/^ (3.5) 

and Ps is a smooth family of pseudodifferential operators with the leading symbol p{x, ^) 
independent of s, so that k, = expHp and p{0,0) — (thanks to (3.2) and p'{0,0) — 
(since k(0, 0) = (0,0)). We shall derive a simple formula for the logarithmic derivative: 
Start with 

hD^ls + PsUjs^O, Uo,s^l, (3.6) 
14 



and recall that Pg and Ut^a commute. Apply hDg to this relation: 

hD^hDsU~J + Ps{hDsU~J = -{hDsPs)U~^, 

which implies 

hD^U^-T,shDsUjs) ^ U^_rjPs + hDj)hDsU~^ ^ -U^_~^{hD,P,)U~^. 
Integrate this from t = to t = t: 

Taking t = 1, we get the promised formula: 

U-\hDsUs) = - f U-t,s{dsPs)Ut,sdt, (3.7) 
Jo 

under the assumption (3.5). 

Let K : neigh (0, R"^) — > neigh (0, R^") be a canonical transformation as in Theorem 
1.3 (so that (1.32) holds), and choose p — po + 0{p^) satisfying (1.33): 

«(p) = expiyp(p) + C»(p~). (3.8) 

Recall that p is uniquely determined modulo 0{p°°) by n and the choice of the quadratic 
form pq with expi^pg = dn{G). 

Let U be an elliptic Fourier integral operator of order associated to the canonical 
transformation k. We look for a pseudodifferential operator P with leading symbol p such 
that 

U = e-'^l^ (3.9) 
Let Pq be a pseudodifferential operator with leading symbol p and put 

C/o = e-*^°/'*. (3.10) 

Let [0, 1] 9 s I— > t/g be a smooth family of Fourier integral operators as above, all associated 
to K (modulo equivalence) and with t/s=o = Uq-, U\ = U ■ We look for a corresponding 
smooth family of pseudodifferential operators Pg, with leading symbol p, such that Ps=o = 
Po, and 

Us = e-'^^'^. (3.11) 

Then P = Pi wiU solve (3.9). 

Since the Us are associated to the same canonical transformation, the logarithmic 
derivative 

Qs = U-^hDsUs, (3.12) 
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will be of order —1 (i.e. 0{h'^^)) with Weyl symbol: 

Q,(p; h) ~ hqs,i{p) + h\sM + ... . (3.13) 

Motivated by (3.7) we shall first look for a smooth family Ps with leading symbol p and 
Ps=Q = -Po) such that 

Qs = - I e'^^^^^{dsPs)e-'^^^'^dt. (3.14) 
Jo 

Denoting the Weyl symbol of Ps by the same letter, 

P,(p; h) = p{p) + hpsM + h'^psM + ... , (3.15) 
we first see that pa,i should solve 

qs,i{p) = - f idsPs,i) o eMtHp)dt + 0{p°^). (3.16) 

JO 

As in the proof of Theorem 1.3, we see that (3.16) has a unique solution dsPs,i (mod 0{p°°)) 
and since po,i is given by the choice of Pq, we get a unique choice of ps,i- 

Proceeding inductively, we assume that we have found pj'"-' with symbol 

m oo 

p^'^Kp; h)^Yl h'PsAp) + E h'PoAp), (3.17) 

j=0 m+1 

where ps,o = P and h^Poj are the terms in the asymptotic expansion of Po(p; h), such that 

- e^^^i"V^(a,pM)e-tP(-)/.^, ^Q^^ h-+'R^+,,,, (3.18) 
Jo 

where Rm+i,s is of order with leading symbol rm+i,s- We just saw how to obtain this 
for m = 1. 

If ^ is a pseudodifferential operator of order 0, we see that e^*^i^VMe-^*^i'"V/^ wiU 
change by an operator of order < — (m + 1) if we modify P^™^ by an operator of order 
< —(m + 1), for instance by passing to P^'^'^^\ It follows that 

To get Pg^™+^^ satisfying (3.18) with m replaced by m + 1, it suffices to have 

r ^^tPi-^/hQ^^p^m+l) _ p{m)y-itPi-^/h^^ = /i^+^P^+l,, + C»(/i-+2) (3.19) 

Jo 

(with the same Rm+i,s as in (3.18)), which gives for the leading symbols 



/ {dsPs 
Jo 



^+l) o exp{tHp)dt = rm+1,8- (3.20) 
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Again this has a unique solution dsPs^m+ii ^ind our induction procedure can be continued 
and gives a solution to (3.14). 

Let Us = e-^^^^l^. Then by construction IJq = Uq, IJ-^hDsUs = U'^hDsUs, and 
Lemma 3.1 implies that Us = Us and in particular that 

U = ^-iP/h^ p=p^ (3 21) 

This gives the existence part of the following 

Theorem 3.2. Let k : neigh (0, R^"^) neigh (0, R^") be a smooth canonical trans- 
formation as in Theorem 1.3 and choose Hj, po as there, so that (1.32) holds. Let 
p E C°°(neigh (0, R^"^); R) be the unique function modO(p°°) of the form p = po + 0{p^) 
with k{p) = expHp{p) + 0{p°^). 

Let U be an elliptic Fourier integral operator of order associated to k. Then there 
exists a pseudodifferential operator P'^{x, hD^', h) with symbol 

P{p;h)^p{p)+hpi{p) + ..., (3.22) 

such that 

U = e-'^'^. (3.23) 

P is uniquely determined modulo "=" and up to an integer multiple of 2'nh by (3.23) and 
the choice of pq. 

It remains to prove the uniqueness modulo "=". Let P'^{x, hD^; h) be another oper- 
ator with the same properties; 

P{p;h)^p{p) + hpi{p) + ..., p = po + 0{p^). (3.24) 

Then we must have k{p) = ex.pH~{p) + 0{p°") and from the uniqueness part of Theorem 
1.3, we conclude that p = p + 0{p°°). 

Put Ps = {1- s)P + sP,0<s< l^so that Pq ^ P, Pi = P and define Us = e'^^-/^. 
For this family, define Qs by (3.12). If P ^ P, let 1 < m < oo be the smallest integer with 

Pm-Pm^O{p'^). (3.25) 

If m = 1, we may also assume that pi — pi is not = to an integer multiple of 2it. From 
(3.14), we see that 

oo 

Qs ^^h^qs,j, 
1 

with qs,j{p) = 0{p°°) for 1 < J < m — 1, and with 

qm,s = - [ {Pm- Pro) o eMtHp)dt + 0{p°°). (3.26) 

^0 
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When m = 1, qm,s is not = to an integer multiple of 27r. ^From the invertibility of the 
map (1.28), we conclude that 

qm := qm,o = qm,s + 0{p^), qm ^ Oip^). (3.26) 

Let Wg, < s < 1 be a smooth family of Fourier integral operators which solves 

Qs = W-^hDsWs, Wo = 1. (3.27) 

If go,s had been rather than just 0{p°^), the Ws would have been pseudodifferential 
operators, so in general they are equivalent to such operators: 

oo 

Ws = R'^ix, hD^; h), where i?,(p; /i) ~ h'rj,s{p), (3.28) 

i=o 

satisfying 
when m = 1 and 

ro,s{p) = 1, rj^s = for 1 < j < m - 2, dsrm-i,s = iqm, 
when m > 2. In other words, 

ro,sip) = e''"^^*^''^ when m = 1, 
Rs{p; h)^l + isgmh"^-^ + C(/i"'), when m > 2. 

In both cases, we have i?i ^ 1, so 

Wi ^ 1. (3.29) 
If we put Us = UqWs, we see that Uq = Uq and that 

Qs = U-^hDsUs = U-^hDsUs. 

By Lemma 3.1 we conclude that Us = Us and in particular, 

UoWi = Ui = Ui. 

Since Wi ^ 1, this contradicts the assumption that Ui = Uq, and the proof of Theorem 
3.2 is complete. # 

Remark. Up to equivalence we have that U is unitary iff P is self-adjoint: 

U*U= l^P* =P. 

Indeed (3.23) gives 

([/*)-! = e-'^*/'', 
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and it suffices to apply the uniqueness statement in Theorem 3.2. 



4. BirkhofF normal forms. 

To get a normal form for the Fourier integral operator in Theorem 3.2, it suffices to 
get the quantized "Birkhoff' normal form of the operator P. For simplicity we shall make 
a non-resonance assumption, and simply recall how this was done in [Sj] in a slightly less 
general setting (in the spirit of works of Bellissard-Vittot, Graffi-Paul and others cited 
there). The extension to the present case is however completely immediate. 

Let P ~ p{p) + hpi{p) + .. be as in (3.22), with p real, p{0) = 0, p'{0) = 0. Put 
Po{p) = |(p"(0)p, p) and let B be the corresponding fundamental matrix, so that 

Po(p) = ^<T(p,Sp), '^B = -B. (4.1) 

Let Hj, — /Xj and possibly be the distinct eigen- values of B. We recall that Theorem 3.2 
was obtained under the assumption (1.32). We add a non-resonance assumption, and for 
that purpose we temporarily change the notation slightly and denote by fij , —fXj , 1 < j < n 
all the eigen- values of B, possibly repeated according to their multiplicity. Assume 

n 

kj/ij = 0, kj e Z ^ ki = .. = kn = 0. (4.2) 

1 

This implies that the /ij are distinct and 7^ 0, so B has the 2n distinct eigen- values fij, —fij, 
^ j which is in agreement with the earlier notation with r = n. 
Notice that (1.32) and (4.2) combine into the single condition 

n 

kj/ij e 27riZ, kj e Z ^ ki = .. = kn = 0, (4.3) 

1 

which does not change if wc modify the choice of the fij by some multiples of 27ri. 

Let ei, .., Cn, /i, ■■, fn G C^" be a basis of eigen-vectors of S, associated to /Ui, .., 
-//I, .., -/In- Then a{ej, e^) = a{fj, fk) = and a{fj, e^) = for j ^ k. We can arrange 
so that 

and then we have a symplectic basis in C^". The corresponding coordinates given 
by C^"^ 9 P = J2i{xjej + ^jfj) will be symplectic, and in these coordinates, we get 

n 
1 

with the Hamilton field 

n 

Hpo = Yl {xj9xj - )• (4-5) 
1 
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If we consider Hp^ : — > V^^^, we see that the monomials x'^^^, \a\ + = m form a 
basis of eigen-vectors and 

Hp,{x^^^)=fi-ia-(3)x'^^^, (4.6) 

where /i = (/Ui, ..,/Un)- The assumption (4.2) impHes that n ■ {a — (3) =0 precisely when 
ct = /3, so if we let the set of resonant polynomials 7^^^^ C Phoj^ be the space of linear 
combinations of all the x"^" = (a;i^i)"i...(x„^^)'^" with 2\a\ — m, we see that Hp^ induces 
a bijection from T't^m/^hom ^^^^ itself. 

We say that u e C°° (neigh (0, R^"^)) is resonant if its Taylor expansion at only 
contains resonant polynomials. Since po is real it is easy to see that the space of resonant 
smooth functions is closed under complex conjugation. We also see that u is resonant iff 
3/ e (neigh (0, C")) with df{p) = 0{p^) such that 

uix) = f{xi^l,..,XnCn)+Oip°°). 

Considering Taylor expansions it is easy to get (cf [Sj]): 

Lemma 4.1 For every v e C°° (neigh (0, R^"^)), 3u G (neigh (0, R^"^)) unique up to a 
resonant function, such that 

Hp^u = v + r, 

where r is resonant. If v = 0{p"^), we can find u,r with the same property. 

As for Hp we only give the corresponding existence statement: 
Lemma 4.2. For every v e (neigh (0, R^'')), 3u e (neigh (0, R^")), such that 

HpU = V + r, 

where r is resonant. If v = 0{p^), we can choose u,r with the same property. 

Notice that since p is real, if v is real, we can take u,r real. The classical Birkhoff 
normal form is then given in 

Proposition 4.3. 3 a smooth canonical transformation k : neigh (0, R^") neigh (0, R^"), 
such that k{p) = p + 0{p^), and 

p o K ^ po + r, 

where r is resonant and 0{p^). 

Proof. If g e C°°(neigh(0,R2'");R), q{p) = 0{p"'+^), with m > 2, then we see that 
expHqip) = p + Hq{p) + 0(p2--i). Let first ^3 = 0(p^) solve Hp{q^) = [p - po) - rs, 
where r3(p) = 0{p^) is resonant. Let K2{p) — Hqaip)- Then 

pMp)) = Pip + HqM + C»(p')) = pM + r3(p) + 0{p^) =: p{p) + rM- 

(We used Lemma 1.1 in [Sj].) Now repeat the argument with p replaced by p and find 
K3 = exp Hq^ et.c. Finally, we choose k with k{p) ~ K2 o K3 o ) o .... See for instance [Sj] 
for more details. # 
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We next review the quantum normal form of a pseudodifferential operator. Let Pq = 
Pq{x, hDrc). If A is (equivalent to) a pseudodifferential operator with symbol a ~ ao{p) + 
hai{p) + we say that A is resonant if every aj is resonant. Since aj is resonant precisely 
when Hpf^aj — 0{p°°) and [-Po, ^] has the symbol ^Hp^a, we see that A is resonant if 
[Po,^] = 0. (Later we shall also recall that A is resonant precisely when it is equivalent 
to a function of the elementary action operators.) 

With p as above, let P = P^ be a pseudodifferential operator with leading symbol p, so 
that P(/0; h) ~ p{p)+hpi{p)+.... Let k be as in Proposition 4.3 and let t/ be a corresponding 
elliptic Fourier integral operator that we choose to be microlocally unitary near 0. Then 
U~^PU = P, where P has the leading symbol p = po + r with r = 0{p^) resonant. 
We drop the tilde and continue the reduction of "P = P" by means of conjugation with 
pseudodifferential operators. We look for a pseudodifferential operator Q = of order 
0, such that 

where R is resonant. Here the left hand side can also be written 

giQpe-^Q = e^^'^op = P + mdgP + i^^^^p + (4.8) 

where the sum is asymptotic in since adgP is of order < —k. We look for Q with symbol 
qo + hqi + .... The leading symbol of iadgP is hHq^p — —hHpqo, so we first choose go so 
that 

Hpqo=pi + ri, (4.9) 

with ri resonant. Then the first two terms in the asymptotic expression of the operator 
(4.8) become resonant. The choice of qi will influence the term in the symbol of 
^iQp^-rQ QYiiy via the term iadgP, and to make the h'^ term resonant, leads to a new 
equation of the same type as (4.9). It is clear that this construction can be iterated and 
we find Q so that (4.7) holds with R resonant. 

If the original symbol P is self-adjoint, then the new P = P = U~^PU" will also be 
self-adjoint and hence have a real-valued symbol. We can then find Q in (4.7) self-adjoint, 
because of the observation that if A, B are self-adjoint, then iad^i? is self-adjoint, so if 
(5, P arc self-adjoint, then all terms of the last expression in (4.8) have the same property. 
Consequently, in each step of the computation, we will encounter an equation of the form 
Hpqk = Pk + fk, with pk real- valued, and we then choose the solution and the resonant 
remainder rk to be real. This means that e"**^ will be unitary. liV = t/e"**^, we finally 
obtain for the original P, that 

V-^PV = Po + R, 

where R is resonant of order with leading symbol r = 0{p^). Summing up we have 

Theorem 4.4. Let p{p) — po{p) + 0{p^) be real-valued with po{p) = ^a{p,Bp), where 
B is symplectically anti- symmetric satisfying the non-resonance condition (4-2). Let P be 
a pseudodifferential operator with leading symbol p. Then there exists an elliptic Fourier 
integral operator V associated to the canonical transformation n in Proposition 4-3, such 
that 

V-^PV = Pq + R (4.10) 
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where R is a resonant pseudodifferential operator of order < with leading symbol = 0{p^). 
If P is self-adjoint, we can choose V to be unitary. 

When applying this to U and P in Theorem 3.2, we notice that 

V-^UV = e-'^"''^/\ (4.11) 

which can be viewed as a quantum normal form for our Fourier integral operator U. 

In the appendix to this section, we review that under the non-resonance assumption 
(4.2), there are real symplectic coordinates xi, .., Xn, Ci, ••, Cn such that 

^hc 

P0{P) = ^{0ij{x2j-1^2j-l + X2ji2j) - /3j(x2j-l6i " X2ji23-l)) (4.12) 
1 

2nhc+''^hr n ^ 

2nhc + l 2nhc+'ihr + l 

where z/j G R are non- vanishing with distinct values of {i/jl, fij > are distinct, and 
aj,Pj > with aj + iPj distinct. We have the corresponding resonant "actions": 

, 1 < J < nhc, (4.13) 



^3=^3^3-: 2nhc + 1 < J < 2nhc + 'T'hr, 

1 

2 



H = o (^j + ^|)' 2nhc + ^hr + 1 < i < 2hc + nhr + ne = n. 



A resonant function is one which can be written /(ti, .., in) + 0{p'^) for some smooth 
function /, and using the simple functional calculus of section 2, we see that a pseudodif- 
ferential operator R of order is resonant iff i? = F(Ii, ..,In', h), where F(t; h) is a classical 
symbol of order and Ij = lJ{x, hD^] h) is the corresponding commuting family of quan- 
tized actions. (We refer to [DiSj] and references there to the original work of B. Helffer 
and D. Robert, for more elaborate functional calculi.) Combining this with Theorem 4.4 
and (4.11), we get 

y-'^UY = Q-'i'F{ii,..,in,\h)/h, (4.14) 

where 

oo 

F(.;/^)^^F,(0/^^ (4.15) 



and 



rihc 2nhc+TJ.hr n- 

^0(0 = Y^^^^^ ~ (^3^n^c+3) + Y ^^^^ + Y ^^^^ + ^^''^)' ^^-^^^ 
1 2nhc+l 2nho+nhr + l 
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Appendix. Review of the real normal form for the quadratic part. 

Here we review some standard facts. See also [Ze2], [It]. Let B : R^" R^"' be 
the symplectically anti-symmetric matrix of Proposition 1.2, case b. We make the non- 
resonance assumption (4.2), so that all the eigen- values of B are simple and 7^ 0. Recall 
from section 1 that they can be grouped into families of 2 or 4 according to the following 

3 cases: 

Case 1. /X > is an eigen-value. Then — /x is also an eigen-value. Let e, / be corresponding 
real eigen- vectors with cr(/, e) = 1, spanning a real symplectic space of dimension 2. A 
point in this subspace can be written p = xe + so that (x,^) become symplectic 
coordinates, and we get po{p) = h{p) = |cr(p, Bp) = ^a{xe + ^f, p,xe — fi^f) = jix^. The 
corresponding resonant action is x^. 

Case 2. p is an eigen-value with Re//, Im^u > 0. Then —p,,]!, —Ji are also eigen- values, and 
we let e, f,e, f be corresponding eigen- vectors. We have 

a{e,e) = 0, a{ej) = 0, a{fj) = 0, (A.l) 



and if e is fixed, we can choose / so that 



a(/,e) = L (A.2) 

e, / and e, / span 2-dimensional complex symplectic subspaces that arc complex conjugate 
to each other, while e, e, /, / span a 4-dimensional symplectic space which is the complex- 
ification of a real symplectic space (of real dimension 4). 
Writing p = ze + (f + we + ujf, we get 

Kp) = ^f^(P> Bp) = iizC + JIwuj. (A.3) 

The resonant action terms are and wu. 
To get the real canonical form, we write 

1 1 

e= ;^(ei +«e2), f = -^{fi - 1/2), 

with Cj, fj real. Using this in (A. 1,2), we get 

a{ej, Cfe) = (T{fj, fk) = 0, a{fj, e^) = 5j^k, (A.4) 

so ei, 62, /i, /2 is a symplectic basis in the real symplectic space mentioned above. 
We also have the inverse relations 

ei = ^(e + e), e2 = ^(e-e), 
/i = ^(/ + 7), /2 = ^(/-7). 
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Write 

2 2 

p = ze + Cf + we + u}J =Y^ XjCj + ^ ^jfj, 

1 1 

so that (x,^) are real symplectic coordinates on our symplectic 4-space. Then 

z = — - 1X2), w = —^{xi + 1x2), 

and using this in (A. 3), we get 

b{p) = a{xi^i + X2^2) - P{xi^2 - X2^i), (A.5) 

with jj, — a + i(3. 

The resonant actions can also be written 

1 

ZC = ^{{Xl^l + X2^2) + i{xi^2 - X2^l)) 
= + ^26) - i{xiC2 - X2Cl))- 

The (resonant) real-valued functions (on the real symplectic 4-space above) which only 
depend on z^,wu> are precisely the functions of xi^i + 0:2^2) xi^2 — X2C1 modulo 0{p°°). 
Notice that these two functions Poisson commute. 

Case 3. /U 7^ is an eigen-value with Re /U = 0. Then fl = —(i is also an eigen- value. If e is 
an eigen-vector corresponding to then e will be an eigen- vector corresponding to /I and 
(T(e,e) G iR \ {0}. Possibly after permuting /i and — and after normalization, we can 
assume that ^ 

Tc(e,e) = 1. 

e,e span a 2-dimensional symplectic subspace which is the complexification of a corre- 
sponding real 2-dimensional space. Let / = ie, so that cr(/, e) = 1. Writing p = ze + C/, 
we see that C, are complex symplectic coordinates, and b{p) = fiz(. 

Write e = -^(^i + ^£2) with ei, 62 real, so that / = -^{€2 + iei) . Then we see that 
ei, 62 is a real symplectic basis. Also notice that 

ei = ;^(e - if), 62 = ;^(/ - ie), 

so if 

p = 2;e -h C/ = a:ei ^62, 
we see that a;, ^ are real symplectic coordinates on our symplectic 2-space and 

b{p) = ^y+e). (A.6) 

24 



The resonant action zQ becomes |(a;^ + times a constant factor. 

5. Parameter dependent case. 

In some applications (for instance when deahng with an energy dependent monodromy 
operator ([SjZw])) our Fourier integral operator will depend smoothly on some real param- 
eter s, and then it may happen that the non-resonance condition is fulfilled for one value 
of s, say for s = but not everywhere in any neighborhood of that point. In this section 
we show that the previous results still apply if we require them to hold only to infinite 
order with respect to s at s = 0. We do this by checking the earlier constructions step by 
step. 

Let A = Ag he real symplectic 2?T,-matrix depending smoothly on s e neigh (0, R), 
such that As satisfies the assumptions of Proposition 1.2, case b. Then log^o can be 
extended to a smooth family of real matrices Bs = log As with e^^ = Ag, '^Bg + Bg = 0. 
(The construction of B = log A can be reformulated by writing B = f{A), where f{z) is 
a suitable holomorphic branch of the logarithm, defined near the spectrum of A. We take 
Bg — f{As) for the same /.) 

Let K*(p), s e neigh (0,R) be a smooth family of canonical transformations with 
K*(0) = and assume that n = kP fulfills the assumptions of Theorem 1.3, so that 

K^{p)^e^pHpo{p) + 0{p^), (5.1) 

where is unique modulo 0{p°") once its quadratic part pg has been fixed in accordance 
with Proposition 1.2.b. We want to extend to a smooth real- valued family p^, with 

K^p) = e^pHps(p) + 0((s,p)~). (5.2) 

Define q' = 0{p^) as in (1.23), so that 

{nydsK' = Hqs, (5.3) 

and consider the problem analogous to (1-24): 

g'(p)= / dgp' oe^ptHp.{p)dt + 0{{s,p)'^). (5.4) 
Jo 

Putting s = 0, we get a unique solution {dgp'^)g=o — 0{p'^) modulo 0(p°°). If we differen- 
tiate k times we get 

[\d^+V) o e^ptHps{p)dt = dU'ip) + Fk{p\ .., a,V,p) + 0((s,p)°°), 

and if .., {d^p^)a=Q = C'(p^) already have been determined, we get {d^'^^p^)g=Q = 0{p'^) 
from this equation. It is then clear that (5.4) has a solution which is unique mod 0((s, p)°°)- 
Let — exp Hps . Then 

{K^rdgK' = {n'TidgK') + 0((s,p)°°), «° = «o. 
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and as in the proof of Theorem 1.3, we see that (5.2) holds. 

We next look at corresponding families of Fourier integral operators and we start 
by extending the equivalence notions of section 2 to the parameter dependent case. If 
Vj C R"^ are open neighborhoods of and Ij C R are open intervals containing 0, we 
say that vj e C°°{Ij x Vj), j = 1,2 are equivalent if they are equivalent in the sense of 
section 2 with Vj there replaced by /-,■ x Vj. Similarly, we define equivalence for symbols 
a^^^ e S^yi^j ^ ^j) ci.nd the corresponding notion for pseudo differential operators. 

Two canonical transformations Kj^s '■ neigh (0, R^"^) neigh (0, R^"^) depending 
smoothly on s G neigh (0, R) with Kj,s(0) = 0, are said to be equivalent, if ki,s{p) = 
K'2,s{p) + C'((s,p)°°). (We write = for the parameter version of equivalence also.) Again 

«l,s = «2,s iff «i,s = K'2,8- 

If we parametrize Kg by a non-degenerate phase (j)s{x,y,9) depending smoothly on s 
(and with (x, y, 9) = (0, 0, 0) corresponding to Ks(0) = 0) then Kg is equivalent to smooth 
family Kg iff we can parametrize Kg by a phase 4>g{x, y, 9) which is equivalent to 4>g{x, y, 9). 

Consider a family Ug = I{ag, (j)g) of elliptic Fourier integral operators as in (2.4), associ- 
ated to a smooth family of canonical transformations Kg as above. Assume 0s (0, 0, 0) = 0. 
We say that ^7^ is equivalent to a second family Ug, if Ug has an associated femily of 
canonical transformations 'Kg with Ks(0) = and we can represent Ug = I{ag,(f>g), with 
0s(O, 0, 0) = and with = a^, (pg = (j)g in the sense of families. (We then have Kg = Kg.) 

Let Ug, Vg be two families of elliptic Fourier integral operators as above, with Uq — Vq 
and Uj^hDgUg = V~^hD8Vs (in the sense of families). Then t/^ = Vg- In fact, let Kg, 'kg 
be the associated canonical transformations and write U~^hDgUs = —Pg, so that Pg is a 
smooth family of pseudodifferential operators of order with real leading symbol Pg{p) = 
C(p2). Then Kg satisfies dgKg{p) = {Kg)^{Hp^{p)). Similarly dsKs{p) = (/5s)*(-f/'~ (p)), 
where pg = pg and kq — kq, so it follows that Kg = Kg. 

Without loss of generality, we may assume that Let Wg be some fixed elliptic 

family of Fourier integral operators associated to Kg, then 

Ug = AgWg, Vg=BgWg, (5.5) 

where Ag, Bg are smooth families of pseudodifferential operators. We get 

U-^hDgUg = W-HAjHhDgAg))Wg + W-'hDgWg, 
Vg-^hDgVg = Wg-\B;\hDgBg))Wg + W-^hDgWg. 

It follows that 

Aj^hDgAg = Bg^hDgBg, Ao = So, 
and then as in the proof of Lemma 3.1, that Ag = Bg and hence that 

Ug = Vg. (5.6) 

We can now prove 

Theorem 5.1. Let Kg : neigh (0, R^"^) neigh (0, R^"^), s G neigh (0, R) be a smooth 
family of canonical transformations with Ks{0) = and assume that kq satisfies the as- 
sumptions of Theorem 1.3. Let U^ = I{as,(f)s) be a corresponding smooth family of elliptic 
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Fourier integral operators of order with (/>s(0,0, 0) = 0. Choose jij and po = Pq as 
there, so that (1.32) holds. Then by that theorem and Theorem 3.2, there exists a real- 
valued smooth function p^ = Pq -{- 0{p^) (uniquely determined mod such that 
K.^{p) — expHpo{p) + 0{p°°) and a corresponding pseudodifferential operator with lead- 
ing symbol p^ , so that = e~'^^ 1^ . (P^ is uniquely determined up to equivalence and an 
integer multiple of 27rh.) 

P° can be extended to a smooth family of pseudodifferential operators so that 

in the sense of families. Moreover, the family P^ is unique up to equivalence for families 
and an integer multipe of2Tvh. The leading symbol Ps satisfies (5.2). 

Proof. Let = {U^)~^hDsU^ be the logarithmic derivative. We first look for solving 

Q' = - f e''^^/^{dsP')e-''^^/''dt. (5.8) 
Jo 

As in section 3, we first determine {dsP^)s=o- Then we can write 

HDsie-'''''/^) = Pt,,e-**^^/'^' (5.9) 

where Rf^s is a well-defined smooth family of 0th order pseudodifferential operators for 
0<t<l,s = 0. We can then differentiate (5.8) once with respect to s and get for s = 0: 

dsQ' = f e^^^^/'^f^i?^,,, a,P*]e-^*^^/'^dt - ^*''''^{dlP')e-'''''/^dt. (5.10) 
Jo " ' Jo 

^From this we determine {&lP^)a=o- Then {dsRt,s) s=o is well-defined in (5.9) and we can 
differentiate (5.10) once more et c and determine {d^P^)s=o for all k. This means that we 
get a solution of (5.8) and we also see that this solution is unique modulo equivalence for 
families. Prom this we also get the uniqueness of in the theorem, for if P* is as in the 
theorem, then it has to satisfy (5.8). 

It remains to show that P* in (5.8) solves (5.7). For that, we put 

ys ^ ^-itP'/h^ 

SO that by (5.8) and the earlier arguments of section 3: 

Qs ^ {V')-\hDsV'). 

Since is also the log-derivative of the family C/* and = V^, we conclude as in (5.6), 
that t/^ = y*, and the proof is complete. 

We end by indicating how to extend the Birkhoff normal form to the parameter 
dependent case. Let 

P^*) ~p(^)(p) + /ipS"^(p) + ..., (s,p) e neigh (0,RxR2"), (5.11) 
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be smooth in (s,p) withp(*)(p) real-valued. Asume that 

pt\p) = \{{p^'^n^)p.p) 

satisfies the non-resonance condition (4.2). Then in suitable complex linear symplectic 
coordinates, we have 

n 

1 

and if we allow the coordinates to depend smoothly on s, we get 

n 

1 

where jif^ depend smoothly on s and jif^ = jij. ^From the appendix of section 4, it 

follows that we can find a real linear canonical transformation Kq*\ depending smoothly 
on s such that 

n 

P^^on^\p) = Y,t^fx,ii, (5.14) 
1 

where the coordinates (a;,^) are now independent of s. After composing P(*) with Kq^ 
we can assume that we have (5.13) with coordintes x, ^ that are independent of s. Notice 
however that the non-resonance condition (4.2) may be violated for s arbitrarily close 
to 0. 

We say that a function r = rs{x, ^) G C°° (neigh (0, R x R^")) is resonant if -H^co) = 

in the sense of families. Notice that this definition does not change, if we replace H (o) by 

H (s) (provided we have (5.13) in s-independent coordinates). Also notice r is resonant iff 

we have = fs{xi^i, ..,Xn^n) for some smooth family/g. The extension of this definition 
to the case of pseudodifferential operators is immediate. We next extend Lemma 4.2: 

Lemma 5.2. For every v = v^^^ e C°° (neigh (0, R^")) there exist u = u^^^ and r = r^*^ 
in C°° (neigh (0, R x R^"')) with r^*) resonant, such that 

Hp^s)U^''^ = v^'^ + r^'l (5.15) 
If V — 0{s^p"^), then we can choose u,r with the same property. 

Proof. Lemma 4.2 gives a solution u^^\r^'^^ for s — 0. Differentiate (5.15) with respect 
to s: 

H^,s)d,u^'^ = - {dj^\u^'^ + dsr^'\ (5.16) 

and put s — 0. Let {dsU^''^^)s=o, {dsr^'^^)s=o be the corresponding solutions to this equation, 
given by Lemma 4.2, then difi^erentiate (5.16) et c. In this way, we get the Taylor series 
expansion of u^^\r^^^ with respect to s, and the lemma follows. # 
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Proposition 5.3. There exists a smooth family of canonical transformations Kg '. 
neigh (0, R^") neigh (0, R^^) with Ks{p) = p + 0{p^) and 

where r^*) is resonant and 0{p^). 

The proof is essentially identical to that of Proposition 4.3. The treatment of the 
operators goes through without any changes, and we get 

Theorem 5.4. Let P^^^ denote also the h-Weyl quantization of the symbol in (5.11). 
Let p^^ be given by (5.13) in the coordinates for which (5.12) holds and let Pq be the 
corresponding quantization. Then there exists a smooth family of elliptic Fourier integral 
operators V = V^^^ associated to k^^^ o (cf (5.14) and Proposition 5.3) such that 

p(s)p^(s) ^ pis) _^ pis) 

in the sense of families, where R^^^ is a resonant pseudodifferential operator of order < 
and with leading symbol = 0{p^). If P^^^ is self-adjoint, then we can choose V^^^ unitary 
(microlocally near 0). 
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